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The partial averaging technique is defined and used in conjunction with the random series im- 
plementation of the Feynman-Kag formula. It enjoys certain properties such as good rates of con- 
vergence and convergence for potentials with coulombic singularities. In this work, I introduce the 
reader to the technique and I analyze the basic mathematical properties of the method. I show that 
the method is convergent for all Kato-class potentials that have finite Gaussian transform. 
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I. INTRODUCTION 

The thermodynamic properties of a monodimensional 
spinless quantum system characterized by the inverse 
temperature /3 = l/ffc^T) are completely determined by 
the canonical partition function 



Z(f3) 



p(x, x\ f3)dx, 



(1) 



where the (unnormalized) density matrix 
p(x,x';/3) = (x\e~ 



>-P H \x') 



can be computed with the help of the Feynman-Kag rep- 
resentation formulaic 



p{x,x';(3) 



E cxp 



-P / V 
Jo 



x (u)+aB° du \ (2) 



for a large class of potentials. In Eq. (0), uiq is the mass 
of the particle, xq(u) is a shorthand for x + (x' — x)u, 
a = (h 2 (3/m ) 1/2 , and 



1 



V2t 



■. exp 



{x-x'f 



2a 2 



is the density matrix of a similar free particle. The 
stochastic element that appears in Eq. {B°, < u < 
1}, is a so-called standard Brownian bridge defined as 
follows: if {B u , u > 0} is a standard Brownian motion 
starting at zero, then the Brownian bridge is the stochas- 
tic process {B u — uBi,0 < u < 1}. Unless otherwise 
stated, in this paper, we shall reserve the symbol E to 
denote the expected value (average value) of a certain 
random variable against the underlying probability mea- 
sure of the Brownian bridge B°. 

The generalization of the Eq. (J2Jl to a (/-dimensional 
system is straightforward. The symbol B® now denotes 
a d-dimensional standard Brownian bridge, which is a 
vector (-B° 1 , B® 2 , . . . ,B° d ) with the components being 
independent standard Brownian bridges. The symbol a 
stands for the vector {a\, oi, . . . , <7d) with components de- 
fined by of = h 2 fl/mo y i- The product aB® is interpreted 
as the (i-dimensional vector of components (TiB®^ Fi- 
nally, x and x' are points in the configuration space R d 
connected by the line xq(u) = x + [x 1 — x)u. 



As emphasized in Ref. ( 3), the success of the Feynman- 
Kag formula in the computation of the thermodynamic 
properties of quantum systems is fortuitously due to an- 
other remarkable result: the Metropolis et alX sampling 
algorithm of arbitrary finite-dimensional probability dis- 
tributions, which lies at the heart of the Monte Carlo 
integration schemes^ This leads to the related but sep- 
arate problem of finding a rapidly convergent sequence 
of finite-dimensional approximations of the stochastic in- 
tegral The main techniques found in literature can 
be classified in two categories: the discrete path integral 
methods and the random series techniques [for reviews, 
see Refs. (@) and (0), respectively]. The latter methods 
lend themselves to various modifications which result in 
convergence for a wider class of potentials V(x) or/and 
improved asymptotic convergence. One such method is 
the partial averaging (PA) technique, which was initially 
introduced by Doll, Coalson, and Freeman2> as a way 
to accelerate the convergence of the "primitive" Fourier 
path integral method (FPI). 9 

As we shall see, the partial averaging method requires 
the Gaussian transform of the potential V for its imple- 
mentation. For real life potentials this is a difficult but 
not impossible task. However, it was generally considered 
that the improvement the technique brings in does not 
warrant the effort of computing the Gaussian transform 
of the potential and therefore, the so-called gradient cor- 
rected partial averaging method was used instead. It has 
been shown that this latter method has general 0(l/n 2 ) 
asymptotic behavior for sufficiently smooth potentials 
and it has been argued that there is not much reason 
to suspect a better convergence rate for the full partial 
averaging method^ However, more accurate numerical 
evidence recently presented in Ref. (0) suggests that the 
full partial averaging method does have in fact a better 
behavior: if the technique is used in conjunction with 
the FPI approach and if the potential is smooth enough, 
the asymptotic convergence is 0(l/n 3 ). The importance 
of the partial averaging method resides also in the fact 
that it acts as a prototypical strategy for improving the 
asymptotic rate of convergence of the random series path 
integral methods. As such, the reweighted random series 
technique^ achieves superior asymptotic convergence by 
simulating the partial averaging approach. 
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In this work, I shall argue for one more property of 
the partial averaging method which is not shared by the 
gradient corrected version and in general by the non- 
averaged methods. More specifically, I shall show that 
the method can be employed for potentials having nega- 
tive coulombic singularities, for which standard discrete 
path integral techniques (and also the primitive random 
series ones) fail to converge iii^ In this respect, it is quite 
surprising that the technique has been scarcely used for 
this purpose, despite the fact that in several instances its 
application for the polaron problemi&i^ and for the com- 
putation of statistical properties of quantum systems^ 
was numerically successful. Though several other meth- 
ods for dealing with the coulombic singularities have been 
proposed (see for instance Ref. Ilfil). I appreciate that the 
advantage of the partial averaging strategy can be best 
emphasized by its ability of handling such systems. 

In this paper, I establish a sufficiently large class of 
potentials for which the partial averaging sequence of ap- 
proximations converges to the correct Feynman-Kac re- 
sult, though I only study the convergence of the density 
matrix and of the partition function. This class includes 
most of the smooth and bounded from below potentials as 
well as most of the potentials having coulombic singular- 
ities. The proofs I perform are of a rather trivial nature, 
as they are direct consequences of well established con- 
vergence theorems from general probability theory. In 
fact, they exploit the martingale property of the partial 
averaging method. Besides their intrinsic value, these 
convergence theorems are important because they set the 
proper mathematical context in which the partial aver- 
aging method should be further discussed or utilized. 



II. FORMULATION OF THE PROBLEM 

A. A chemically relevant class of (scalar) potentials 

A sufficiently large class of potentials for which the 
Feynman-Kac, formula J2Jl and its multidimensional ana- 
logues hold is the so-called Kato class, which we define 
below. If 



Ivl d = l, 
g{y) = { Hh\\- X ) d = 2, 
\\y\\ 2 - d d>3, 



then the Kato class Kd is made up of all measurable 
functions / : R d -> K such that 



lim sup 

Q J-° xEM d J \\x-y\\<a 



\f(y)g(x- y)\dy = 0. 



(3) 



We also say that / is in K l d oc if 1 of 6 Kd for all bounded 
domains D C M. d . We say that the potential V(x) is of 
Kato class if its negative part V- = max{0, — V} is in 
Kd, while its positive part V+ — maxjO, V} is in K d oc . 
In these conditions, as shown in Refs. fl^) and |lg|) . the 
Feynman-Kac formula holds. Moreover, Th. B.7.5 



of fRef. Il8|) shows that the density matrix p{x,x';0) is 
continuous on R d x R d x (0,oo), while Th. B.6.7 of the 
same reference shows that for a given (3 > 0, the density 
matrix is uniformly bounded in the variables (x, x'). 

A remarkable theorem due to Aizenman and Simon^S 
gives an alternative definition for the Kato class Kd- 
More precisely, Theorem 4.5 of Ref. (|19f) says that V S 
Kd if and only if 



lim sup E 



\V(x + aB u )\du 



= 0, 



(4) 



where E denotes the expectation value against the d- 
dimensional Brownian motion B u . Inverting the order 
of integration in Eq. (JIJ and remembering that B u is a 
Gaussian distributed variable of variance u, we obtain 
the equivalent condition 

lim sup f du f (2Tiu)- d/2 e-^ 2/{2u) \V{x+az)\dz = 0. 

e i° x£M d JO JK d 

(5) 

We leave it for the reader to perform the substitutions 
u' = u/e and then z' — z/y/e in successive order and 
prove the following reformulation of the condition given 
by Eq. ©: 



du / (2wu) 



-d/2 e -\\zr/(2v 



\V(x + ayflz)\ dz = 0. 



lim sup e 

x£R d JO 

(6) 

In the Appendix, we shall use Eq. in the proof of 
Th. Eland Eq. © in the proof of Th. 0] respectively. 

As far as the chemical physicist is concerned, the Kato 
class of potentials is sufficiently general. It contains for 
instance the coulombic potential as it appears in elec- 
tronic structure calculations. For another example, the 
ab initio intermolecular potential computed at the level of 
the Born-Oppenheimer approximation cannot have sin- 
gularities worse than the coulombic ones and therefore 
it is of Kato class. However, we do not consider certain 
empirical potentials which are not of Kato class, as for 
example the Leonard- Jones potential. Nevertheless, this 
can be brought into the Kato class if the unphysical r~ 12 
singularity is removed by truncation or by other approx- 
imations. 

Let us anticipate a little and also demand that the 
potential V have finite Gaussian transform. The reader 
may read ahead in the next subsection and see that this 
condition is natural for the proper definition of the partial 
averaging method. More precisely, we require that 



\V\ a (x) 




° XP ( ~5I^2 ] \ V ( X + Z )\ <°° 



(7) 



for all xeff'andaeR'f. [In this paper, R + = (0, oo).] 
Certain properties of the potentials having finite Gaus- 
sian transform are given by Theorem [3] of the Appendix. 
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From the thermodynamic point of view, only the diag- 
onal density matrix p(x, x; 0) is of interest. Moreover, in 
order to have a physically relevant statistics, the condi- 
tion 



< Z{0) < oo; V/3 > 



(8) 



must hold, but the inequality JSJ is not a requirement for 
the results obtained in this paper to be valid. In practice, 
the condition is achieved by the addition of a constrain- 
ing potential, which is usually a continuous and bounded 
from below function (thus still in the Kato class). The 
constraining potential is intended to simulate, for exam- 
ple, the effect of the container in which a reaction takes 
place. A sufficient condition for the quantum partition 
function to be finite is that the analog classical parti- 
tion function be finite. This follows from the following 
inequality: 



Proposition 1 Set 



-0V(x) 



dx. 



Then, Z({3) < Z c i{(3). 

Proof. By Jensen's inequality and Tonclli theorem, 



Z(f3) = 



1 



v / 2na 2 



< 



V27R7 2 



dxE due~^ x+ ° B °^ 



= — ^=E Cdu f dxe-^ x+ " B ^ = Z d ((3). □ 
V2na 2 Jo Jm 

As stated, Proposition |Q remains true for multidi- 
mensional systems. In this paper, we shall perform the 
proofs only for monodimensional systems. The reader 
should notice that our arguments are purely measure the- 
oretic, in fact irrespective of the dimensionality of the 
physical systems. On the other hand, in the chemical 
physics literature it is customary to perform the analy- 
sis in "monodimensional" notation. I consider that the 
mathematician will have little trouble generalizing the 
proofs, yet the chemist may find it hard to accommodate 
a more complicated notation. 



B. The partial averaging strategy 

In this section, we shall give a short review of the par- 
tial averaging method for monodimensional systems (the 
multidimensional generalization is straightforward). For 
a more complete discussion, the reader should consult 
Ref. (0). The most general series representation of the 
Brownian bridge is given by the Ito-Nisio theorem^ the 
statement of which is reproduced below. Assume given 
{Afc(r)}fc>i a system of functions on the interval [0,1] 
which, together with the constant function Aq(t) = 1, 



makes up an orthonormal basis in L 2 [0, 1]. Let A k (u) 
denote the primitives 



Afe(u) = / A fc (r)dr 
Jo 

of the functions X k (u). If £1 is the space of infinite se- 
quences a = (di, a2, . . •) and 



P[a] = U /i(o fc ) 



(9) 



fe=i 



is the (unique by the Kolmogorov extension theorem) 
probability measure on Q such that the coordinate maps 
a — ► a k are independent identically distributed (i.i.d.) 
variables with distribution probability 



Li(a k G A) 



1 



e~ z / 2 dz, 



then 



B° u (a) = ^a fe A fc (u), 0<u<l 



(10) 



(11) 



k=l 



is equal in distribution to a standard Brownian bridge. 
Moreover, the convergence of the above series is almost 
surely uniform on the interval [0, 1]. 

Using the Ito-Nisio representation of the Brownian 
bridge, the Feynman-Kac formula © takes the form 



P{X / X '\ P) = [ dPa;exp(- 3/V 
Pfp{x,x';(3) J n I J 



a^a k A k { 



k=l 



An 



(12) 



The independence of the coordinates a k , which physically 
amounts to choosing those representations in which the 
kinetic energy operator is diagonal, is the key to the use 
of the partial averaging method. Denoting by E„ the av- 
erage over the coefficients beyond the rank n, the partial 
averaging formula reads 

P p n A {x,x'-p) = 
p fp (x,x';[3) 

x exp <^ - (3 E r , 



dju(ai) . . . / d/i(o n ) 



V 



x (u) + o'^2a k A k (u) dtt>.(13) 



k=l 



Assuming that the Fubini theorem holds (this is proved in 
the next section) , one may invert the order of integration 
in the exponent and compute 



l 

2„ / V[x {u) + aBl(a)]du 
lo 
fi 

E n V[x (u) + aB°{a)]du 

o 

/ V Utn [x (u) +a y2<i k A k (u)]du, 

J fe=l 



(14) 
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where 

Vu,n(y) 



: exp 



The function T^(u) is defined by 

oo 



2ra(u) 



V(y + z)dz. 

(15) 



fe=n+l 



i(l-u)-^Mu) 2 



fc=l 



(16) 



[Again, the reader is refered to Ref. (.7) for additional 
explanations. It is customary to use the notation r^(u) 
to mean the square of T n (u).] To summarize, we define 
the n-th order partial averaging approximation to the 
diagonal density matrix by the formula 



= / dju(ai) . . . / dp,{a n ) 



Pf p (x,x';(3) 



xex.pl -(3 / V Uin x (u) + (j^2a k A k (u) du>. (17) 



III. 



k=l 



MARTINGALE PROPERTY AND 
CONVERGENCE RESULTS 



This section establishes the martingale property of the 
partial averaging method. One may notice that we ob- 
tain some important convergence results without actually 
doing much work other than citing some well-established 
theorems. The chemical physicist will probably be more 
interested in the Corollary 1, which is for that matter 
presented separately. 

On the set f2 of sequences a = (0,1,0,2, ■■ ■), consider 
the tr-algebra generated by the finite-dimensional Borel 
sets Too = cr(U„>o^ r n), where T n = a(a 1 ,a 2 , . . . , o„) and 
J-q = {0,J7}. By construction, {T„} n >o is a filtration. 
Also, the probability measure dP[a] introduced in the 
previous section is, of course, defined over Too so that 
the default probability space we refer to in this work is 
(SI, Too, P). If / is an integrable random variable on SI, 
we shall sometimes denote the conditional expectation 
E [/ \T n ] simply by E„/. 

To continue with the introduction of the notations, we 
define 



U n (x, x', (3; a) = / V u ,; 

Jo 



•?'o 



(u) + a^a k A k (u) 



k=l 



du 



and 



U 00 (x,x',(3;a) = / V 
Jo 



x (u) + a^a k A k (i 



fc=i 



du. 



The variables x, x' , and (3 are interpreted here as pa- 
rameters and, just as a reminder, we separate them by a 



semicolon from the "true" variable a. By construction, 
U n (x, x' , f3;a) is T n measurable, while Uoo{x, x' , (3; a) is 
Too measurable. Let us prove that U n (x,x',(3;a) — 
E n [Uoo(x, x' , (3; a)] for P-almost every a. As shown in 
the previous section, this boils down to proving 



E, 



V[x (u) + aB°(a)]du 



E n V[x (u) + aB°{a)}du 



P-almost surely. The relation follows from the Fubini- 
Tonelli theorem provided that 



E„ \ V[x (u) + aB°(a)}\du 



(18) 



is finite P-a.s. Since the above integrand is nonnegative, 
it is enough to show that its P-expectation is finite. Let 



1 



-.e- z ' 2 /^dz 



denote the respective Gaussian measure on K. Using the 
Tonelli theorem to invert the order of integration, one 
computes 



y.\ / E n \V[x (u) + aB°(a)}\du 



E\V[x (u)+aB°(a)]\du (19) 
du / \V[x a (u) + z]|d^r o(u) (z) < 00. 



The last expression is finite for all {x,x';j3) € M 2 x K + 
by Th. ^ of the Appendix. Eq. l[T9")l also shows that 
U 0o (x,x',f3;a) is P-integrable for all (x,x';(3). Then, 
standard theorems from martingale theory show [see 
Theorem 5.7 of Ref. ©]: 

Theorem 1 For all (x,x',(3) e M 2 x R + , the sequence 
U n {x, x' , (3; a) is a martingale adapted to the filtration T n 
and is a.s. and L 1 convergent to E [Uoo(x, x' , (3; a)\T oo\ = 
Uoo(x, x' , (3; a). 

Let us define 

X n (x, x', (3; a) = pf p (x, x';(3) exp[-^ U n (x, x' , (3; a)] 
and 

Xoo(x, x', (3; a) = pf p (x, x'; (3) exp[-^ Uoo(x, x' , (3; a)]. 
Then, we have 

p(x,x';/3)=E[X o0 (x,x',L3;a,)] (20) 

and 

pF A (x,x';(3)=E[X n (x,x',0;a)], (21) 
respectively. 
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Theorem 2 (PA Convergence Theorem) For all 

(x, x',/3) € R 2 x R + , </ie sequence X n {x, x' , (3; a) is a 
submartingale adapted to the filtration T n and is a.s. 
and L 1 convergent to X co (x 1 x' , (3; a). 

Proof. We notice that X n (x,x' , (3;a) is the exponen- 
tial of a martingale. Thus, by the conditional Jensen's 
inequality [see page 225 of Ref. lf2l|) ]. we have 



X n <E[X n+1 \T n ] KElX^n] 



(22) 



The above inequality establishes the submartingale prop- 
erty because, as mentioned in Section II. A, 

p(x,x';l3)=E[E[X 0o \F n }] ^EX^ 

is uniformly bounded in the variables (x, x') for all (3 > 0. 
An elementary proof of this assertion is given in the Ap- 
pendix (see Th.[5J). Finally, the a.s. convergence follows 
directly from Th. ^ while the L 1 convergence follows 
from the Dominated convergence theorem and the in- 
equality □ 

We define the n-th order partial averaging partition 
function by the formula 



Using the symbol f to mean "monotonically increasing 
to," a direct consequence of Th.[3]is the following. 

Corollary 1 As n — > oo, 

p? A (x, x'; (3) T p(x, x 1 ; f3) and Z? A ({]) | Z((3). (23) 

Proof. The pointwise monotonic convergence of the den- 
sity matrix is a direct consequence of the submartingale 
property and of the L 1 (fi, P) convergence of the partial 
averaging method. Then, the convergence of the par- 
tition functions follows from the Monotone convergence 
theorem. □ 

IV. SUMMARY 

In this paper, I presented the basic properties of 
the partial averaging method. I demonstrated that the 
method can be employed for a quite general class of po- 
tentials by proving several convergence results of interest 
for the chemical physicist. In particular, I proved that 
the PA method is convergent for most of the potentials 
having negative coulombic singularities. The value of the 
convergence theorems deduced in the present paper con- 
sists of the fact that they establish the mathematical con- 
text in which the partial averaging method should be uti- 
lized or discussed. I also anticipate that the martingale 
property will play an important role in establishing the 
asymptotic rates of convergence for different partial av- 
eraging schemes and, in fact, it may explain the superior 
asymptotic behavior of these methods. 
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APPENDIX A 

The following theorem consists of well known facts 
about the Gaussian transform and I present it here for 
ease of reference. 

Theorem 3 Let f : M. d — ► R be a Borel measurable 
function, let a = (ao, . . . , ay) € R^J_, and let \\a\\ — 
maxKKd cti. Consider the d- dimensional Gaussian mea- 
sure 



id let 



d/i Q (z) = Yl 



F(x,a) 



1 



(Al) 



\f(x + z)\d/i a (z). 



be defined on R d x 

such that F(y,rj) < oo and let D = I 
Then the following are true: 

a) f is locally integrable. 

b) F(x, a) < oo for all (x, a) G D. 
c) 



Assume there is (y, 77) e R d x Rl 

t d x n?=i(o,»*). 



G(x, a) 



fix + z)dfi a (z) 



is well defined, continuous and infinitely differen- 
tiable in both arguments on D. 

d) limi^ii^g G{x, a) = f(x) a.e. More strongly, 



lim 

llall^O 



\f{x + z) - f(x)\dp, a (z) = a.e. 



Theorem 4 Let x and x' be arbitrary points in R d , let 
(3 > 0, and let B® be a standard d-dimensional Brownian 
bridge on [0, 1]. Pick some arbitrary a = {a\, . . . , ad) € 
R^ and set y = x' — x. If V G K l d oc and V has finite 
Gaussian transform, then 



[ E\V(x + yu + crB u )\du 
Jo 



du \V(x + yu + z)\ d/i ro („)(z) < 00, (A2) 



where Tq(u) = u(l — u)(af, . . . , Og) and the Gaussian 
measure dp, a (z) is defined by the relation fAlp . 
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Observation. We stated this theorem separately be- 
cause its proof depends upon the dimensionality of the 
problem. More precisely, if the system is monodimen- 
sional, one may use Tq (it) < a 2 to show that the integral 
121) is smaller than 



du- 



1 



/o - u) 

By Theorem 01c), the integral 



\V{x + yu + z)\du a {z). (A3) 



\V(x + yu + z)\ d^icr(z) 



as a function of u is continuous on [0, 1], thus bounded. 
Then, by the integrability of [u(l — it)] -1 / 2 , it follows 
that the integral (|A3|I is finite. However, this reasoning 
is not valid for higher dimensions because [it(l — u)]~ d / 2 
is not integrable for d > 2 and the additional condition 
V G K d oc is needed. 

Proof of the theorem. The equality i|A2fl was discussed 
in the text. From Th. |3|c), it follows that 



/ 



\V(x + yu + z)\du ro(u) (z) 



as a function of u is continuous on all compact intervals 
[e, 1 — e] with < e < 1/2, thus bounded and integrable. 
It is then enough to show that 



Ie(x,y) = du \V(x + yu + z)\dfi ro{u) (z) 



du \V(x + z)\ dnr (u)(z - yu) < oo 



for all x and y and small enough e. This is so because 
the integral over the end [1 — e, 1] can be shown to equal 
I e (x + y, —y) by the change of variable v! = 1 — u. We 
shall prove the above inequality in two steps. 
Step 1. In the first step, we prove the inequality 



I e {x,y) < 2 d exp 



4a 2 



/ du 
lo 



\V{x + z)\d^^{z). (A4) 
The inequality < u < 1/2 implies 



as a quadratic function of Zi is attained at Zi = 2yiU and 
has the value —yfu 2 . Therefore, 

(zi ~ yiu) 2 _ z 2 /2 + zf/2 - 2z l y i u + yfu 2 



2o 2 u 



> 



2a(u 

9 2 9 9 

zf _y±}t > _jt_ y[_ 

Aa 2 u 2o 2 ~ A<j 2 u 4a. 2 ' 



where we used again the condition u < 1/2. Replacing 
the last inequality in Eq. (| A5|> . we obtain 



n 



L = \ [ y/2lT<J 2 u(l - U) 



exp 



(zi - yiuf 



2a 2 u{\ - u) 



< 



i=l 



: exp 



4a 2 u 



and, consequently, the inequality given by Eq. (|A4I) is 
proven. 

Step 2. By the results from the first step, it suffices to 
show that the last integral in Eq. (|A4|> is finite. By ap- 
propriate transformation of coordinates, we may rewrite 
this last integral as 



du / \V(x + z)\dn aV ^(z) 



du / \V(x + z)\du aV ^(z) 



e I du 
10 



V(x + zaV2e) du^(z) 



Even more, the integral over M. d in the last relation can 
be restricted to the ball \\z\\ < 1. Indeed if ||z|| > 1, we 
have 



(2ttm) 



-d/2 



exp 



1*11' 



= u d / 2 exp 



1 

1 

u 



Notice that 



2u J I 2 

x (2^)- d/2 exp(-||z|| 2 /2) < u-^e^e" 1 '^) 
x (2^r d/2 exp(-||z||72). 

,1/2, -d/2 , 



M = supe i/2 u- d/2 exp[-l/(2u)] < oc 

ji>0 

Then, a little calculus and Th. |3]c) show that 



n 



i 



l J[ [ yj2-KCF 2 u(l - U) 



exp 



jzj - y t u) 2 
2a 2 u(l - it) 



< 



2 d / 2 ]J 



■. exp 



(zi - yiuf 



2a fu 



i= i 1 V^a 

On the other hand, the minimum of the expression 
z 2 /2 - 2ziy { u + y 2 u 2 



(A5) 



e I du 

>0 J\\z\\>l 



< eM Q 



V{x + zaV2e) d^(z) 
V(x + zaV2e) d/ii(z) < oo. 



To conclude the theorem, we only need to prove that 



e j du 

'o J \\z\\<i 



V(x + zaV2e) d/i^(z) < 
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for e small enough. Pick an arbitrary r\ > 0. Remem- 
bering that e < 1/2 and taking x such that ||x|| < r), we 
notice that in order to compute the integral over z in the 
above formula, we only need to know the potential over 
the ball D of radius r\ + maxKi^Cj centered at origin. 
Therefore, if we set Vd — IdV, then 



e / du 

'0 J \\z\\<l 



< e 



V(x + zav2e) dp,^(z) 
f du f V D {x + zcrV2e) dn^{z) (A6) 



for all x such that ||a;| < r\. Since D is bounded and 
V e K\ oc , it follows that V D € K d and then Eq. © 
guaranties that there is eo > such that the last integral 
in Eq. I|A6J| is uniformly bounded for all x. Consequently, 



e / du 



\V{x- 



z\\<l 



/ 2eo)| dfi^{z) 



is bounded for all x such that ||x|| < n. Since rj is arbi- 
trary, we are done. □ 

Theorem 5 Assume V is a Kato-class potential. Then 
there is Mp > a constant depending upon the inverse 
temperature /3 such that p(x, x'; (3) < Mp for all (x, x') € 
W l x W l . 

Observation. The proof of this theorem will show why 
the Kato class is the natural class for the treatment of 
Feynman-Kag semigroups. Most of the arguments used 
in the proof are borrowed from Aizenman and Simon. 19 

Proof of the theorem. If V— denotes the negative part 
of V, we notice that 

Eexp j-/^ V[x (u)+o-B° u ]\ 

<Ecxpj/?y VL[a;o(u)+(TB°]| 

so, without loss of generality, we may assume that V is 
of class Kd- The proof of this theorem is organized in 
three steps, each step reducing the problem to a simpler 
statement. 

Step 1. In the first step, we prove that it suffices to 
show that 



sup 



p(x, x'; f3)dx' < 00 



(A7) 



for all (3 > and V £ K d . 

In this part of the proof, it is convenient to denote the 
density matrix by pv(x,x';f3), the index V indicating 
the potential from which the density matrix is derived. 
For the proof, we need two well-known properties of the 
density matrix py(x, x'\0): it is symmetrical 

p v (x,x';f3) = p v (x',x;/3) 



and it satisfies the Trotter product rule 

p v (x,x';/3) = / p v (x,y;P/2)p v (y,x';P/2)dy. 



These two properties can be established by direct com- 
putation starting with the definition of the Brownian 
bridge. The first one is a consequence of the symme- 
try of the standard Brownian bridge, that is {B®_ u : < 
u < 1} is also a Brownian bridge and is equal in distri- 
bution to {B° : < u < 1}. The Trotter product rule 
is a consequence of the Markov property of the Brown- 
ian motion B u entering the definition of the Brownian 
bridge. The simple proofs are left to the reader. 

Now, the Cauchy-Schwartz inequality gives the esti- 
mate 



p v (x,x';P) < 



p v (x,y;(3/2) 2 dy 



Pv (y,x';f3/2) 2 dy 



1/2 



1/2 



Taking the supremum over x and x' and using the sym- 
metry of the density matrix, one concludes that 

supp v (x,x';P) < sup / pv{x,y;l3/2) 2 dy. (A8) 

x.x' x JK d 

Again by the Cauchy-Schwartz inequality, 

py(x,x';/3/2) 2 



E cxp i — 



V 



x {u) + —B° u 



p fp (x,x'; [3/2)1 
2 

du 



< Eexpl - 



2V 



zo(«) + 



du 



Pfp (x,x';[3/2) 2 < (n^=?) 



PV {x,x>-i3/2) 2 < (n^=H 



Next, we combine the last equation with the bound 

p fp (x,x';/3/2), 

to obtain the inequality 

p 2V (x,x';[3/2). (A!)) 
Substituting Eq. (|A9|) in Eq. (|A8|) . one obtains 
suppy(x,x';/3) < ( TT j=z ) sup / p 2 v(x,y;f3/2)dy 

x,x> V 71 "^ / x "'K' 1 

and the claim of Step 1 is concluded because (3/2 > 
and 2V € K d . 

Step 2. Simple transformations of coordinates show 
that 

p(x,x';(3)dx' = f dp 1 {z)Ee-' 3 tt> v \ x+zu ' J+ ' jB > u 
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and from the very definition of the Brownian bridge, we 
learn that zu + B® is in fact a Brownian motion B u start- 
ing at zero. Thus, 



p(x,x';(3)dx' = Ee-^» v ( I+ffJJ ") du . 



For the remainder of the proof, E stands for the expecta- 
tion value with respect to the underlying measure of the 
standard Brownian motion B u . 

In this second step, we use the Markov property of the 
Brownian motion to show that if there is eo > such 
that the inequality 



sup 



< oo 



holds for all e < eo, then Eq. (|A7|) also holds. 

Let 9 and r be some positive real numbers such that 
6 + t = 1 . We break the integrand in the above equation 
in two parts 

e -{SQ V(x+oB u )du _ e -/3j° V(x+vB u )du 
xe -Pl e 1 V(x+aB u )du _ e -/}tfV(x+aB u )d U 
xe -P!;V{x+aB e+u )du_ 

Using the Markov property, we learn that the expected 
value of the above integrand conditioned on the random 
variables B u with u £ [0, 0] is 

-P /„" V(x+aB u )du^t e -p JJ" V(x+aB e +aB' u )du 

Here, the symbol E' denotes the expected value against 
the standard Brownian motion starting at zero B' u , Brow- 
nian motion that is independent from B u . The above 
conditional expectation is smaller or equal than 

e -pf Q e V(x+aB u )du sup E e -/3/ T V(x+<rB u )du 

x 

where the prime sign becomes superfluous and is there- 
fore dropped. Taking the total expectation and then the 
supremum over x, we obtain the inequality 

supEe-^"^^")^ < su P Ee-^» v ^ nB ») du 

X X 

X 

A simple inductive argument then shows that 



supine 



PfiV(x+*B u )du < < J sup ]E e -/3/o 1/ "l / (^B u )d« 



(A10) 

for all n > 1. Clearly, the claim of Step 2 is concluded 
because the right-hand side of Eq. IjAlOf) is finite for all 
n such that 1/n < 6q. 

Step 3. In this final step, we prove that there is eo 
small enough such that 



sup 



E e -^/o V{x+uB u )du 



for all e < eo. Eq. Q allows us to pick some eo > such 
that 



supE 



(3 I \V\ (x + aB u )du 
o 



< 1/2 



(All) 



for all e < eo. Next, we consider the inequality 

supEe-^o v(x+„B u )du < supEe pJ*\v\(x+*B u )du 

X X 

00 oh, r p£ ~| x 

= su P ^ttE / \V\(x + aB u )du <^^ fe ,(A12) 
x k=i Uo 



k=0 



where 



Ak = sup 



= /3 fc E f 

jQ<s x <— 



\V\ (x + aB u )du 



\V\(x + aB Sl ) 

<s k <e 

V\ (x + aB Sk )dsi . . . ds k . 



(A13) 



Notice that the term by term integration of the first series 
appearing in Eq. (|A12(I is guaranteed by the Monotone 
convergence theorem. The last equality in Eq. (|A13(I fol- 
lows by symmetry arguments. 

To construct a bound for the terms , we first condi- 
tion on the random variables B u with u € [0, Sk-i] and 
use the Markov property of B u to show that this condi- 
tional expectation has the value 

H k ' x [ \V\(x + aB Sl )...\V\(x + aB Sk _ 1 ) 

J0<S 1 <...<S k _ 1 <£ 



E'/3 



\V\ (x + aB Sk _ 1+ aB' u )du 



o 



dsi...ds fe _i(A14) 



Eq. i|All|) shows that the quantity in the square brackets 
is bounded by 1/2. Therefore, the conditional expecta- 
tion given by Eq. (| A 14(1 is bounded by 



-J 



fe-i 



\V\(x + aB Sl ) 

0<si<...<s^_i<e 

|V| (x + crB Sk _ 1 ) dsi . . .dsfe_i. 



(A15) 



Taking the total expectation in the Eqs. (| A 14(1 and (|A15|) 
and then the supremum over x, we learn that 

A k < A k -i/2, 

from which the inequality Ak < l/2 fc follows by induc- 
tion. Substituting this last inequality in Eq. <|A12|) . wc 
obtain 

oo 

supEe^/o ^+- B «) d " < J- ^ = 2 < oo 



fc=0 



< oo 



for all e < eo and the proof of Step 3 and of the theorem 
is concluded. □ 
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